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The objective of this study is to assess the impact of a dense-phase treatment on the hydrodynamic 
description of granular, binary mixtures relative to a previous dilute-phase treatment. Two theories 
were considered for this purpose. The first, proposed by Garzo and Dufty (GD) [Phys. Fluids 14, 
146 (2002)], is based on the Boltzmann equation which does not incorporate finite- volume effects, 
thereby limiting its use to dilute flows. The second, proposed by Garzo, Hrenya and Dufty (GHD) 
[Phys. Rev. E 76, 31303 and 031304 (2007)], is derived from the Enskog equation which does account 
for finite- volume effects; accordingly this theory can be applied to moderately dense systems as well. 
To demonstrate the significance of the dense-phase treatment relative to its dilute counterpart, the 
ratio of dense (GHD) to dilute (GD) predictions of all relevant transport coefficients and equations 
of state are plotted over a range of physical parameters (volume fraction, coefficients of restitution, 
material density ratio, diameter ratio, and mixture composition). These plots reveal the deviation 
between the two treatments, which can become quite large (>100%) even at moderate values of the 
physical parameters. Such information will be useful when choosing which theory is most applicable 
to a given situation, since the dilute theory offers relative simplicity and the dense theory offers 
improved accuracy. It is also important to note that several corrections to original GHD expressions 
are presented here in the form of a complete, self-contained set of relevant equations. 



I. INTRODUCTION 

Polydisperse, rapid solids flows are quite prevalent in both nature (i.e., landslides, avalanches) and industry (i.e., 
pharmaceutical processing, high- velocity fluidized beds), though much remains to be understood. Perhaps most 
importantly, due to differences in size and/or material density of each particle species, polydisperse mixtures are well- 
known to exhibit particle segregation, also known as de-mixing [1-4]. Such behavior has no monodisperse counterpart. 
Thus, continuum models developed for monodisperse flows cannot be used to predict the segregation of unlike particles 
which occurs in polydisperse systems. Consequently, an accurate continuum model of a polydisperse solids mixture 
lends itself to a variety of non-trivial applications, such as the design of coal gasificrs for energy production. 

The scope of the current study pertains to binary mixtures of inelastic grains (negligible fluid phase) engaging in 
instantaneous, binary collisions (rapid flows). Numerous previous contributors have proposed continuum theories for 
such systems (for recent review, see Ref. Q), and the application of these theories has led to a better understanding 
of the mechanisms by which de- mixing occurs pit TgI— Tl6l] - Nevertheless, the improvement of existing models remains 
an active area of research due to differences in the derivation process. More specifically, one or both of the following 
simplifications have been incorporated in the vast majority of previous models: (i) Maxwellian velocity distribution 
[l7l - l2l| . and/or (ii) an equipartition of energy [22|, [23[. (The theories proposed by Rahaman et al. [2(| as well as Iddir 
and Arastoopour [2l| assumed a Maxwellian velocity distribution between unlike particles only.) The aforementioned 
assumptions are strictly true for systems of perfectly elastic spheres in a uniform steady state [24|, but not so for 
^ ■ inelastic grains. Furthermore, numerous studies have demonstrated the influence of non-equipartition on species 
segregation ll2l-[T5l]. Two continuum models have been proposed over the past decade in which neither of the 
above conditions is assumed. The first theory, developed by Garzo and Dufty (GD) [25T - I27T ] is based on the Boltzmann 
equation, and is thus applicable to dilute flows only. The second theory, developed recently by Garzo, Hrenya and 
Dufty (GHD) [28|, [29[, instead uses the Enskog equation as its starting point, making this theory applicable to 
moderately dense systems as well. Hereafter, the acronyms GD will be used to refer to the former, and GHD will be 
used to refer to latter. Each theory gives rise to a set of zcroth-order closures known as equations of state, as well as 
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constitutive relations for first-order contributions to fluxes, or more specifically the associated transport coefficients. 
The equations of state and transport coefficients are functions of the hydrodynamic variables: number densities (rij), 
mass-based mixture velocity (U), and number-based mixture granular temperature (T). Although the predictions for 
the equations of state and transport coefficients from the two theories are expected to match at low volume fractions, 
a non-negligible difference is expected at higher concentrations, though the level of discrepancy between the two has 
not yet been reported for polydisperse systems. 

To build on the previous contributions, the focus of this work is to analyze binary mixtures, where the two particle 
species differ in mass and/or size. Motivation for this study is threefold: (i) to assess importance of dense-phase 
corrections to hydrodynamic description of mixtures proposed by Garzo and co-workers [28], [29| compared to the 
previous dilute-phase description [25|, and more specifically to determine rules-of-thumb for the volume fraction at 
which such dense-phase descriptions become non- negligible, (ii) to examine the behavior of the GHD equations of 
state and transport coefficients over a range of physical parameters, and (iii) to provide a complete, self-contained set 
of the GHD expressions, including several corrections (see the Appendix [AJ for the expressions given in the original 
GHD contribution [28|, [2^] . This latter goal also provides an opportunity to display the expressions in a form more 
suitable for computational purposes. 

To accomplish the first two objectives, the equations of state and transport coefficients were evaluated over a range 
of volume fractions, coefficients of restitution, and mixture properties (diameter ratio, mass ratio, and volume fraction 
ratio). The results indicate that the discrepancy between transport coefficients and equations of state predicted by 
each theory at a volume fraction of <f> = 0.1 can vary from a factor of f .05 to a factor of 10. As the volume fraction 
becomes fairly dense ((f) = 0.5), the predicted discrepancy increases to a factor of at least 1.7 and as large as a 
factor of 120. Hence, though the derivation of the constitutive relations for a dilute flow and the resulting constitutive 
expressions are simpler than its moderately dense counterpart, the difference between the two theories is non-negligible 
at low to moderate volume fractions. In the upcoming sections, a complete, self-contained set of the GHD constitutive 
relations for the mass flux, heat flux, pressure tensor, and cooling rate are given in Sections |TIJ Also, a quantitative 
comparison between the GHD and GD predictions for the transport coefficients and equations of state illustrates stark 
differences between the dilute and dense treatments (Sections IIIII and IIV[) . The paper is closed in Sec. [V]with a brief 
summary of the main results obtained here. 



II. ENSKOG KINETIC THEORY FOR MASS, MOMENTUM AND HEAT FLUXES AND EQUATIONS 

OF STATE OF A GRANULAR BINARY MIXTURE 

The mass, momentum, and granular energy balances for the GHD theory for an s-component mixture are given in 
Table 1, along with the corresponding flux laws. Each balance equation is expressed in terms of the hydrodynamic 
variables (rii, U, T), along with the following quantities: cooling rate (£), mass flux (ji), heat flux (q), and pressure 
tensor (P). Constitutive expressions for these latter quantities, also given in Table 1, are in terms of (zeroth-order 
cooling rate), ( u (transport coefficient associated with first-order cooling rate), Djj (mutual diffusion coefficients), Df 
(thermal diffusion coefficients), Df, (mass mobility), A (thermal conductivity coefficient), D q: i (Dufour coefficients), 
Lij (thermal mobility), p (pressure), r\ (shear viscosity), and k (bulk viscosity). Also, refers to the external force on 
a particle of species i. To fully close the set of equations, these quantities must be cast in terms of the hydrodynamic 
variables. The equations needed to obtain closures for each expression are detailed in the following subsections; 
corresponding equation numbers are also listed in Table 1. 



A. Mass and heat fluxes 



We consider a binary mixture (s = 2) of inelastic^ smooth, hard disks (d = 2) or spheres (d = 3) of masses mi 
and 77i2, and diameters o\ and oi. The inelasticity of collision among all pairs is characterized by three independent 
constant coefficients of normal restitution an, «22, and u\i = ol%i, where is the coefficient of restitution for 
collisions between particles of species i and j. For moderate densities, it is assumed that the velocity distribution 
functions of each species are accurately described by the coupled set of inelastic Enskog kinetic equations [10, HJ ■ 
This set of equations has been recently solved in Refs. [2|| [29| by means of the Chapman- Enskog method [24| and the 
constitutive equations for the mass ji and heat q fluxes have been obtained up to the Navicr-Stokcs order (first-order 
in the spatial gradients). In the absence of external forces (F^ = 0), the forms of ji and q are given, respectively, by 

m?ni miTO2?i9 t 

ji = — £>iiVlnni " " £>i2Vlnn 2 -p£>fVlnr, 2.1 

P P 



q = -T 2 D 9i iVlnni - T 2 D ?i2 Vlnn 2 - AVT, 



(2.2) 
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where p = mini + 7712x12 and 7ij refers to the number density of species i. While the diffusion coefficients Dij and the 
thermal diffusion coefficient Dj have only kinetic contributions, the transport coefficients D q i and A associated with 
the heat flux have also collisional transfer contributions. Expressions for these transport coefficients in terms of the 
coefficients of restitution, the parameters of the mixture (masses, sizes and composition), and concentration (solid 
volume fraction) have been obtained in Ref. [29j by using the leading terms in a Sonine polynomial expansion. 

Table 1: Hydrodynamic description of a granular mixture from GHD theory. 



Balance Equations 



Mass 


Drii „ TT 1 „ . 
— + ?i,V-U = V-ji 

Dt rrii 


Momentum 


p-jx =-v-p + X] n » F » 

i=l 


Granular 
Energy 


n = Vq P : VU nTC+ T> V-b+> J 

i=l 1 i=l 1 



Flux Laws 



Mass 



n ,m , — I ), ,Y In rij — pDfV In T 



Eq. 2.1 



Heat 



q = - AVT - T 2 D^V In n, - ^ ^ L y F, 

i=l i=l j=l 



Eq. 2.2 



Pressure 



P a p = pS a/ 3 - V (V Q t//3 + VpU a - | V • US a0 ) - kV • UScp Eq. 2.22 



Cooling Rate 



C - C (0) + C«v • u 



Equations of State 


Transport Coefficient 


Zeroth-Order 
Cooling Rate 


£(0) 


Eq. 2.38 


First-Order 
Cooling Rate 


Cu 


Eqs. 2.40, 2.41, 2.42 


Pressure 


P 


Eqs. 2.24, 2.25 


Shear Viscosity 


V 


Eqs. 2.27, 2.28, 2.30, 2.35 




Bulk Viscosity 


K 


Eqs. 2.34, 2.36 


Mutual Diffusion 


Dij 


Eqs. 2.3, 2.7, 2.8 


Thermal Diffusion 


DJ 


Eqs. 2.3, 2.6 


Thermal conductivity 


A 


Eqs. 2.4, 2.13, 2.16, 2.17 


Dufour 




Eqs. 2.4, 2.13, 2.15, 2.20 


Mass Mobility 


D F - 




Thermal Mobility 


Lij 



The transport coefficients depend upon the temperature, concentration, and composition of the mixture, as well 
as the masses, diameters, and coefficients of restitution. To present the expressions of these coefficients in a compact 
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form, it is convenient to consider their dimensionless forms: 

n * _ m\v Q mim 2 i/ Tif _ pv Q T 

V u ~ ~0T~ L>11 > u i2- — ~jjf< — ^ 12 ' v i ~ ^j 7 1 ' ^ 6 > 

D*.= 2 K+^j 2 (m 1 + m 2 H 

«'* d + 2 n 9 ' + 2 nT v ; 

Here, i^o = n<7 i2 lv o is an effective collision frequency, <7i 2 = (01 + er 2 )/2, and uo = y/2T/m is a thermal velocity where 
m = (mi + m 2 )/2. Thus, the results given throughout the remainder of the paper arc given in terms of the mass 
ratio m 1 /m 2 , the size ratio <J\/ct2, the species number fraction Xi = rii/n, the concentration ncr 2 , and the coefficients 
of restitution an, a 22 , and a±2- 



Al. Mass flux transport coefficients 



In a binary mixture, since ji = — j 2 , the mass flux contains three relevant transport coefficients: Dn, D\2, and Df. 
The remaining coefficients D 22 , -D 2 i, and D% are given by the identities 



D 2 i = -D u , D22 = --D12, Dl = -Df. 

m 2 m 2 

The expressions of the reduced coefficients Df*, -D ll5 and D\ 2 can be written as [29| 

Dj* = (vh-CY 1 {xili- ^ + ^^rX 1 na 2 l [x lX ii(<J 1 /<J2) d li(l + an) 



+2a; 2 Xi2(o'i2/o'2)' i Mi272(l +0:12)] |, 



(2.5) 



(2.6) 



C* 



Df* <9C (0) mi dp 971 
"1-5 ^ n i"5 h 7i + 

+ /rfx xin^ X! Xn( cr «/ cr 2) d Mi'i(l + an) 

rfi V 2 / f=i 

2dif + ri£ — 1 lui 

oni ri\ 





mi \ 




7f 




m e J 



mi c*7f 

n ^"5 — 

m^ ani 



(2.7) 



D 



12 



£>T* ^(0) mi 9p g 7l 
n 2-^ ^ni- r-TH — 

n d/2 2 

-fJd\ Xina 2 ^2xu(cr u /a2) d M ei (l + a u ) 



dT 

1=1 
- ( fni 



I621 + ni—r 1 1 U 2 

on 2 n 2 



mi d-fi 
— ng- — 

me cm 2 



(2.8) 



In these equations, pi = raiiii, 7; = Tj/T, Q* = Q^" 1 /i/q, p* = pj (nT), Xij is the pair distribution function at contact, 
Mij = rrii I (rrii + rrij ) , and 



V2V d -D/ 2 



D 



71 72 



-T- — Xi2(l + ai 2 ) — h — — (X1M12 + X2M21) 

* M12 A'hi 1 



1/2 



(2.9) 



The partial temperatures T\ and T 2 are determined from the condition ([^ = = C^°\ where the expression of 
is given by Eq. (|2.38l) . Moreover, an explicit form for Xij f° r disks (d = 2) and spheres (d — 3) is given in the 
Appendix [B] 
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The parameters Iuj are chosen to recover the results derived by Lopez dc Haro et al. for elastic mixtures |32| . 
These quantities are the origin of the primary difference between the standard Enskog theory and the revised version 
for clastic collisions [33|. They are zero if i = £, but otherwise are not zero. They are defined through the relation [29[ 



2 

e=i 



d In Xi 

dflj 



9 'fij ' 



D, 



- 2njXij(Ti 



(2.10) 



where B 2 = -K d l 2 jdY (d/2) [T refers to Gamma function, such that B 2 = ^ for d = 2 (disks) and B 2 = 2n/3 for d = 3 
(spheres)] and Hi is the chemical potential of species i. Taking into account Eq. (|2.10|) . the nonzero parameters im 
and 7i22 appearing in Eqs. (2.11) and (2.12) are given by 



h 



21 



TB 2 n 2 af 2 xi2 



»i 



3hx 
dni 



-T 



T,n 2 



- 2 



"lofxii n \ a \ gXn "i ^Xi2 
"2^12X12 n-2crf 2 xi2 9rii X12 dni ' 



(2.11) 



-^122 — 



TB 2 cr( 2 xi2 \dn 2/ T 



affix dxii n 2 dxi2 



vf 2 Xi2 dn 2 X12 dn 2 



(2.12) 



Explicit forms of Hi f° r disks (d = 2) and spheres (d = 3) are given in the Appendix [51 

A2. Heat flux transport coefficients 

The heat flux requires going up to the second Sonine approximation. Its constitutive equation is given by Eq. (|2.2[) 
where the transport coefficients D q j and A have kinetic and collisional contributions 



ii,. iK., • x = x k + x c 



(2.13) 



The corresponding reduced forms D k \, D c * i , X k * , and A c * are defined as 

n fc,c* _ 2 (mi + m 2 )v ^ c fc _ 2 (mi + m 2 > t 



(2.14) 



The kinetic parts D k \ and X k * can be written, respectively, as 

72 



D k \ = d*n +d*, + ( - 

,,1 9,21 ^ Mi2 



D k q * 2 =d* q22 +d* ql2 



7i 72 



M u M 2 i 



X 2 L>12, 



(2.15) 



A fc * = AJ + A 2 



7i 



12 



12 



21 



u l i 



(2.16) 



where the expressions for the (dimensionless) coefficients d* 4 • and A* are displayed in the Appendix [B] In Eqs. (|2.15[) 
and (|2.16p . the coefficients Df* , and D\ 2 are given by Eqs. (|2.6|) . (|2.7|) . and (|2.8p . respectively (first Sonine 
approximation) . 

Let us consider now their collisional transfer contributions. In the case of the thermal conductivity, A c * is given by 

M 



T d/2 



2 2 



2d(d + 2)r(f) 

+ {mx+m 2 )Dj* 



XE 1 *^/' 72 )^^* 1 + "«) {[(5 ~ a ijWij - (! - a ii) M 3i\ *j 



i=l j=l 



((5 - on^Ma - (1 - ay)-^) + ^ ((3 + (Xij)Mji - (7 + c^)M y -) 

7TZ 7 7TZ a 



16 W -(trxa/traj^ay/aiajq 



S-^/tt mi + m 2 



(2.17) 
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where 



C*- = (Oi + flj)" 1/a (^^)" 3/2 { 2 4 + + (*< + ft) [(ft + OjWijMji + /%(1 + Mjj)] } 



+i (l - ay)^ - M«) 



7 » ~ 
ftft 



3/2 



Here, ft = (m.;T '/ 'mTj) and /3y = Myft — Mji$i. In the case of the coefficients £>£* , they can be written as 

= + % . ^"2 = + D c * 22 , 
where the coefficients D°* j have the explicit forms 

3 tt<V2 2 



(2.18) 



(2.19) 



2d(d + 2)r(f) 
+(mi + m 2 )xjD* j 



d 2 J2 x P {cj lp /a 12 ) d Xlp M lp {l + a jp ){ [(5 - ct y )Af. ip - (1 - 0^)-^*] 



p=i 
7 P 



((5 - a ip )M ip - (1 - ai P )M pj ) 



+— ((3 + a ip )M pi - (7 + a ip )M ip ) 



32 



where 



(ft + P r 1/2 ( W~ 3/2 k>/Mft + e p ) l -e t e p [1 + Mei^+M ^ 



7r mi + m 2 W J 

51n7 p 



(2.20) 



<91nn 7 - 



3/2 



3 ft In 7 P 
2ft + ft,dlnn,; 



(2.21) 



B. Pressure tensor 



The overall constitutive relation for the pressure tensor is a combination of the zeroth (PW) and first-order (P^ 1 -*) 
contributions, which is given by 



where 



p = p(0) + p(!) 



(2.22) 



(2.23) 



The zeroth-order contribution to the pressure tensor is proportional to the mixture granular pressure p. The equation 
of state that defines p is given by 



p = p k + p c 

where the kinetic (p k ) and the collisional (p c ) contributions are [28| 

n d/2 2 2 

p k = nT , p c = ^ jjr p k nai ^ ^ x i x j {ai i /a 2 ) d M j i (1 + ay) Xijli- 

V 2 J i=l j=X 

The constitutive equation for the pressure tensor , proportional to the velocity gradients, is 



(2.24) 



(2.25) 



(2.26) 



Here, r\ is the shear viscosity and k is the bulk viscosity. The coefficient r\ has kinetic and collisional contributions 
while k only has a collisional contribution k c (and so, vanishes for dilute gases) 



T] = rf + rf, k = k c . 



(2.27) 
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The kinetic part n k is 
where the partial contributions iff can be written as 



k nT k* 

m = —Vi- 
vo 



The reduced coefficients 77^* are given by 



2(2r 2 * 2 - 20m - 4Tf 2 r? 2 



C* 2 - 2C*(r* 1 + r| 2 ) + 4(^7-2*2 - Tf^) ' 

2(2^ - 2C*)?? 2 - 4^^ 
C 2 _ 2C(r * i + T * 2 ) + 4{t * iT * 2 _ T * 2Tk y 



where the expressions of the (reduced) collision frequencies r 4 * can be found in the Appendix A of Ref. [2£ 
(|2.30[) . we have introduced the quantities 



Vl = x lJl + E n + E 12, V2 = x 2l2 + E-22 + E 2 1, 



where 



r d/2 



Eij = 



d{d + 2)T 



7jrncr%XiXj((Jij /a 2 ) d XijMji(l + ay) 



MjiQonj - 1) + —-,) ~ 4M ii ( 7< - 7i) 

The collisional contribution rf to the shear viscosity and the the bulk viscosity n have the forms 

nT „ 



where 



c " T c* 

77 = rj , ft 

^0 



2 2 



wo 



2 ?r d/2 ^ 

F7I) d(d + 2) ™ CT 2 £ J2 X i (°* J fa j Mji (1 + tty ) Jfc* 



d + 2 



ft , 



(2.28) 
(2.29) 

(2.30) 

(2.31) 
In Eq. 

(2.32) 



(2.33) 



(2.34) 



(2.35) 



d 2 T ( mi + m 2 



2J ^ a^i^j 
<=i 3=1 



m» + mo (j. 



2(/ 



Xy(l +a, 3 ) 



1/2 



(2.36) 



It must be remarked that the predictions of the shear viscosity 77 compare quite well with Monte Carlo simulations 
of a heated granular binary mixture, even for strong dissipation [34( . 



C. Cooling rate 

The overall cooling rate can be written as the sum of the zeroth-order (C^ -*) and first-order contributions (( u ) 

C = C (0) + C«V • U. (2.37) 

The zeroth-order cooling rate of each species (C^) defines the rate of kinetic energy loss for that species, and is given 
by the following relation 



C (0) = C (0) 



47r (d-l)/2 



-V 



1/2 



(1 + ay) 



Ma , 

l--fO- + <*ii)- L a- L 



(2.38) 
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As shown in Eq. (|2.38|) . the zcroth-ordcr cooling rate for each species is equivalent (i.e., £i°' ) = C^)- Because Eq. 
(j2.38[) is an implicit expression that depends on individual species granular temperatures, the following equation is 
needed 

nT = mTi +n 2 T 2 . (2.39) 

The expressions given in Eq. (|2.38[) and Eq. (|2.39[) form a set of two non-linear algebraic equations that can be solved 
for 0\ and 9 2 (using the relation 0, = to;T '/mT,-), and then species temperatures T\ and T 2 can subsequently be found. 
The equation of state defining was first proposed by Garzo and Dufty [35| . 

At first order in gradients, there is a contribution to the cooling rate from V • U. The proportionality coefficient Q u 
is a new transport coefficient for granular fluids. Two different contributions can be identified 

C« = C (1 ' 0) +C (M) . (2.40) 

The coefficient t^ 1 ' ) is given by 

Q_rf/2 2 2 

C( = "^fTiT"^ EE^ M ^>2) c Wl - a%H. (2.41) 
The contribution (V' 1 ) can be written as 

C(M) = ^fTiT S S ^^ r ^K-M 2 ) d - 1 X«M i ,(l - a^.)^ 3/2 ^ 1/2 (^ + ^)- 1/2 e* D , (2.42) 

where 



2(2V 2 * 2 - 3C*)ei,D - 4^ 12 e 2 , 



9C* 2 - 6C*(^n + ^22) + 4(^*^2*2 - ^2*1) ' 
* = g(2^i - 3C*)e 2 , g - 4^ 2 *iei ,p 

-2,d 9C 2 _ 6Cirii + r22) + 4(^2*2 - r 12 r 2 i) ' 

3 collision frequencies ?/>,*■ have 

given by 



(2.43) 



(2.44) 

Here, the collision frequencies ip^ have been determined in the Appendix A of Ref. [29| and the coefficients e^u are 



n d/2 2 

+4(13 + 2d + <.)n,,}.\l,, - 1n.\/;,W, ill, + 6j)(l + ay) 2 

+15Mf i ej 2 (6 i + 9^(1 + a t] f] . (2.45) 

The results displayed along this section give the explicit forms for the equations of state, the transport coefficients 
and the cooling rate of a moderately dense granular binary mixture. The corresponding expressions for a low-density 
binary mixture can be easily obtained from their dense forms by taking the limit na^ — > 0. These explicit ex pre ssions 
are displayed in the Appendix [C] and agree with those previously derived from the Boltzmann equation [25l |27|| . 



III. QUANTITATIVE APPROACH: COMPARISON OF DILUTE AND DENSE-PHASE EXPRESSIONS 

FOR HARD SPHERES 

In order to assess the importance of dense-phase corrections to the continuum theory for rapid granular flows 
of binary mixtures, the equations of state and transport coefficients obtained from the GHD and GD theories were 
compared over a range of volume fractions and coefficients of restitution for a given set of mixture properties (diameter 
ratio, size ratio, and volume fraction ratio). To illustrate the differences in a straightforward manner, each quantity 
is examined as a ratio of the GHD value (dilute through moderately dense) to the GD value (dilute limit), giving 
rise to a non-dimensional quantity. These non-dimensional ratios were plotted as functions of volume fraction and 
coefficients of restitution, holding all other mixture properties constant. Representing the transport coefficients and 
equations of state in this manner reveals the relative magnitudes of the dense- and dilute-phase predictions. Recall 
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the complete set of equations of state and transport coefficients for GHD theory are given in Table 1 Cu, Dij, 

D i> D ij> \ D q,U L ij> V, V, «)■ 

It is important to note that some transport coefficients (Lij, Dfj) were not considered in the dilute theory (GD), 
and thus these quantities are not considered here. Moreover, two of the transport coefficients, namely Cu and k, are 
zero in the dilute limit, and thus the corresponding ratios of the moderately dense (GHD) value to the dilute (GD) 
value diverge. Accordingly, only the GHD predictions of these quantities are shown. Thus, the comparison between 
the GHD and GD theory predictions presented here involves the seven remaining quantities: C , P, Vi Dij, Df ', A, 
D q ,i- 

A. Mixture parameters 

The continuum description of a binary mixture of inelastic hard spheres (d = 3) is a function of the following 
dimensional parameters: species masses (mi, 1712), species diameters (tri, 02), species I volume fraction (4>\), overall 
volume fraction (0 = (j)i+4>2), and coefficients of restitution (an, 022, 0^12 = ^21)- (Note that the number densities and 
volume fractions are related by 4n = niiraf/6.) The subscripts 1 and 2 denote the two species in the binary mixture. 
For purposes of simplicity, the coefficients of restitution have been assumed to be the same for all combinations of 
collisions (i.e., an = «22 = ^12 = a). In terms of the ratio of moderately dense (GHD) to dilute (GD) predictions for 
each quantity, the parameter space is reduced to the following dimensionless inputs: mass ratio (mi/m2), diameter 
ratio {oxjo-i), overall volume fraction (0), volume fraction ratio of species I relative to the total (<f>i/4>), and coefficient 
of restitution a. Hereafter, the ratio 4>i/4> will be referred to as the (mixture) composition of species I. Recall that 
the GHD and GD theories allow for a non-equipartition of energy, and thus several of the aforementioned closures 
(see, for example, Eq. (|2.f 5[0 involve the species granular temperatures, T\ and T2. It is important to point out that 
these quantities are not hydrodynamic variables (i.e., they do not require the solution of species energy balances; for 
a detailed explanation, see Ref. [29j ) and instead are determined by the set of equations defining the zeroth-order 
cooling rate (Eqs. (|2~5T| and (pOS)) ). 

B. Parameter space evaluated 

Table 2 summarizes the three cases (equal size and different mass, equal mass and different size, and different size 
and mass) used to compare the GHD and GD theories, and the wide ranges of input parameters used in each case 
study. Though the transport coefficients and equations of state may vary quantitatively from case to case, the general 
trends show little variation. For sake of brevity, the upcoming section will focus on one representative case, namely 
that of different size and equal material densities (i.e., different mass) in order to quantify how the newly acquired 
GHD predictions differ from the dilute-phase counterpart (GD). 



Table 2: Range of input parameters used in analysis of binary mixture via GHD theory 



Case Description 


0"l/<7 2 


mi/m2 


<f>l/<f> 


a 


Equal Sizes 


1 


1-10 


0.25-0.75 


0.50-0.99 


Equal Masses 


1-10 


1 


0.25-0.75 


0.50-0.99 


Different diameters, different masses 


2 


0.10-10 


0.50 


0.75 



C. Case presented: Different-sized particles with equal material densities 

Many industrial and natural granular flows are comprised of one material (i.e., same material density), but different- 
sized particles. In the case presented here, the diameter of species 1 was twice that of species 2 (i.e., a\jai = 2), and 
both species had the same material density (i.e., m\jrai = 8). For the sake of consistency, the composition of each 
species was held constant at 50% by volume for this analysis (i.e., <fii/<f> = 4>2/4> = 0.5). The ratio of GHD to GD 
predictions of each quantity evaluated was plotted over a range of volume fractions from dilute to moderately dense 
((f) = 10 _8 -0.5) while holding the coefficient of restitution constant. Also, each quantity was varied over a range of 
coefficients of restitution from relatively inelastic to nearly elastic (0.5-0.99) while holding the overall volume fraction 
constant. The results of this case study are presented in the upcoming section. 
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IV. RESULTS AND DISCUSSION 

The overall objective was to analyze each transport coefficient and equation of state over a range of parameters for 
the newly-developed GHD theory. By comparing these results to the predictions from the dilute (GD) theory, it was 
possible to demonstrate the need for a moderately dense-phase correction, as detailed below. 








FIG. 1: Zeroth-order cooling rate: ratio of moderately dense (GHD) to dilute (GD) predictions as a function of (a) overall 
volume fraction and (b) coefficient of restitution. 



0.1 0.2 0.3 0.4 0.5 0.6 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 




FIG. 2: Transport coefficient associated with first-order cooling rate: moderately dense (GHD) predictions as a function of 
(a) overall volume fraction and (b) coefficient of restitution. See legends presented in Figure 1. 



A. Cooling rate: zeroth-order and first-order contributions 

As indicated by Figure 1, the dense-to-dilute ratio of the zeroth-order cooling rate (C^) is much more sensitive 
to changes in volume fraction than coefficient of restitution. Such behavior can be explained via the dependency of 
the zeroth-order cooling rate (Eq. I2.35| on the pair correlation function at contact, \ij 

(Eqs. (jBTif and (|BT7)) V This 

factor accounts for the volume exclusion effects between like particles (xn) and unlike particles (xi2)- in the dilute 
limit, the spatial correlation factor equals one (i.e., xn = X12 = !)• When the zeroth-order cooling rate of GHD 
theory is then divided by its dilute counterpart, the resulting function is strongly dependent on the spatial correlation 
factor. Because Xij is sensitive to changes in overall volume fraction, it is then reasonable that the dimensionlcss 
zeroth-order cooling rate ratio exhibits the same sensitivity. More specifically, the results shown in Figure la indicate 
that C (0) predicted by GHD theory is more than 5 times greater than its dilute counterpart for a fairly dense system 
(4> = 0.5) and more than 2 times greater for <fi = 0.3. Even at = 0.2, a discrepancy of 27% is found between the 
dilute- and dense-phase predictions. 

Unlike the zeroth-order contribution to the cooling rate, the transport coefficient associated with the first-order 
contribution is zero in the dilute limit. Therefore, a ratio comparison of the dense-to-dilute predictions is not possible. 
The results given in Figure 2 represent the first-order contribution to the cooling rate (which is non-dimensional), 
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which approaches zero as volume fraction diminishes. As evident from this figure, £ u is quite sensitive to changes in 
both volume fraction (Fig. 2a) and coefficient of restitution (Fig. 2b). Also, the results for this case indicate that the 
magnitude of the first-order contribution increases as the system becomes denser and less elastic. 




FIG. 3: Pressure: ratio of moderately dense (GHD) to dilute (GD) predictions as a function of (a) overall volume fraction and 
(b) coefficient of restitution. See legends presented in Figure 1. 




a 



FIG. 4: Shear viscosity: ratio of moderately dense (GHD) to dilute (GD) predictions as a function of (a) overall volume 
fraction and (b) coefficient of restitution. See legends presented in Figure 1. 



B. Momentum flux: pressure, shear viscosity and bulk viscosity 

Now moving on to results associated with the momentum flux, Figure 3 indicates that the dense-to-dilute ratio of 
granular pressure is more sensitive to changes in volume fraction (Fig. 3a) than coefficient of restitution (Fig. 3b). 
Also, this ratio increases monotonically with both volume fraction and coefficient of restitution. For a moderately 
dense system (</> ~ 0.4), GHD theory predicts that the granular pressure is about 5 times greater than dilute (GD) 
theory. Even at lower volume fractions (4> ~ 0.1), the moderately dense-phase prediction is greater than its dilute 
counterpart by 40%. 

Shear viscosity, results of which are given in Figure 4, behaves in a similar manner to granular pressure (Fig. 3). 
A monotonic increase is exhibited with respect to both volume fraction (Fig. 4a) and coefficient of restitution (Fig. 
4b). The GHD prediction is about 5 times larger than the GD prediction for moderately dense systems (0 ~ 0.4). 
However, the discrepancy at lower volume fractions (<f> ~ 0.1) decreases to approximately 5% (Fig. 4a). 

As mentioned previously, the bulk viscosity is zero in the dilute limit. Therefore the results for bulk viscosity, 
given in Figure 5, are those obtained from the moderately dense theory (GHD) alone, instead ratios of dcnsc-to-dilute 
predictions. Furthermore, these GHD-bascd bulk viscosities are non-dimensionalized according to Eqs. (|2.34[) and 
(|2.36[) . It is evident from this figure that the prediction of bulk viscosity via GHD theory increases significantly in 
magnitude as the system becomes moderately dense, whereas little variation results from changes in particle elasticity. 
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a 

FIG. 5: Bulk viscosity (non-dimensional): moderately dense (GHD) predictions as a function of (a) overall volume fraction and 
(b) restitution coefficient. See legends presented in Figure 1. The dimensionless inputs are as follows: mi/m2 = 8, ai/<J2 = 2, 
and 4>i/4> = 0.5. 




FIG. 6: Mutual diffusion (-Dn): ratio of moderately dense (GHD) to dilute (GD) predictions as a function of (a) overall 
volume fraction and (b) coefficient of restitution. See legends presented in Figure 1. 



C. Mass flux: Mutual diffusion, thermal diffusion 

The mutual and thermal diffusion coefficients (D^j, Df) are elements of the constitutive equation for the mass flux. 
Based on the identities given by Eq. (|2.5j) . the dimensionless mutual diffusion can be described by two quantities 
(Dn/Dn^iiute and D22/ D22, dilute), whereas the dimensionless thermal diffusion can be described by a single quantity 
(Df /Df di i utc ). Note that the dilute (GD) mutual and thermal diffusion coefficients presented by Garzo and Dufty 
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FIG. 7: Mutual diffusion (-D22): ratio of moderately dense (GHD) to dilute (GD) predictions as a function of (a) overall 
volume fraction and (b) coefficient of restitution. See legends presented in Figure 1. 
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[25j are defined using different spatial gradients than those used in the dense (GHD) theory and shown in Eq. (12.11) . 
Nonetheless, a conversion is made such that both dense and dilute theories use the same representations for the fluxes, 
namely those shown in Scc.|n]of this paper, thereby ensuring an apples-to-apples comparison. 

An examination of the dense-to-dilute ratio of the mutual diffusion coefficient elements (Figures 6 and 7) reveals a 
more complicated behavior of these quantities. Because -Dn,diiute approaches zero for inelastic systems (at a ~ 0.52), 
the coefficient of restitution was varied between 0.6 and 0.99 (Fig. 6b). As shown in Figure 6a, Du/Dn. dilute is non- 
monotonic with respect to the volume fraction in less elastic systems (a = 0.5), reaching a maximum ratio between 
the dense and dilute predictions of 20 at a volume fraction of 0.37. As the system becomes more elastic, Dn/Dn t dilute 
shifts from positive to negative. A change in the sign, as well as magnitude, between dense and dilute pred ictions of 
the mutual diffusion coefficient Z? n may provide insight into counter intuitive species segregation fi\, M. \l4. Il5l|. 

The results for -D22/-D22, dilute, which are displayed in Figure 7, reveal increasing discrepancies between predictions 
as volume fraction increases and restitution coefficient decreases. In other words, GHD and GD theories display a 
larger discrepancy in denser, less elastic systems. For a relatively inelastic and dense system (a = 0.5 and <j> = 0.5), 
the GHD prediction is about half of its dilute counterpart (Fig. 7a). However, it is significant to note the minor 
differences that exist between the dense and dilute predictions for the mutual diffusion coefficient D22 near the elastic 
limit (a = 0.9) over a range of volume fractions from <\> = 10~ 8 to 0.5 (Fig. 7a, -D22/-D22, dilute ~ !)■ GHD and 
GD theories display a larger discrepancy in denser, less elastic systems. For a relatively inelastic and dense system 
(q = 0.5, <fi = 0.4), the moderately dense-phase theory prediction is about half of its dilute counterpart (Fig. 7a). 
Comparing dense- and dilute-phase predictions for the individual elements D\\ and -D22 shows the relative importance 
of each contribution to the mutual diffusion. At a moderately low volume fraction and high coefficient of restitution 
(<p = 0.1, a = 0.9), the discrepancies for GHD and GD theory predictions are about 70% and 5% for Dn and D22, 
respectively. The dilute theory does not consider the finite size of the particles, which is the main difference between 
dense and dilute predictions. The discrepancy between Dn and -Dii,diiute is larger than the discrepancy between D22 
and D22, dilute because Dn is directly related to the size of species 1, whereas -D22 is proportional to the size of species 
2 (recall a 1 /a 2 = 2 for the case examined). Neither dilute quantity contains species size, therefore, the self-diffusion 
coefficient of a relatively large particle compared to its dilute counterpart will be greater than that of its smaller 
counterpart. 




FIG. 8: Thermal diffusion: ratio of moderately dense (GHD) to dilute (GD) predictions as a function of (a) overall volume 
fraction and (b) coefficient of restitution. See legends presented in Figure 1. 



The results for thermal diffusion (Figure 8) indicate that the ratio of dense-to-dilute predictions is extremely 
sensitive to changes in volume fraction compared to the coefficient of restitution. These general trends were also 
observed in the cooling rate and momentum flux relations (Figures 1, 3-5). The quantity /D^ diiutc is nearly linear 
when plotted as a function of volume fraction, regardless the restitution coefficient (Fig. 8a). In a moderately dense 
system (0 = 0.4), the results of Fig. 8a indicate that the dilute (GD) theory prediction of Df is 5 times larger than 
predicted by GHD theory. At a much lower volume fraction of 0.1, the dilute (GD) theory prediction is larger than 
the moderately dense- phase (GHD) theory prediction by 20% (Fig. 8a). 



D. Heat flux: thermal conductivity, Dufour coefficients 



Heat flux is characterized by the thermal conductivity A and the Dufour coefficients D q ^. Figure 9 shows that 
the dense-to-dilute ratio of thermal conductivity increases monotonically with respect to both volume fraction and 
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FIG. 9: Thermal Conductivity: ratio of moderately dense (GHD) to dilute (GD) predictions as a function of (a) overall volume 
fraction and (b) coefficient of restitution. See legends presented in Figure 1. 






FIG. 10: Fig. 10. Dufour coefficient (D qi i): (a) ratio of moderately dense ( GHD) to dilute (GD) predictions as a function 
of overall volume fraction (b) dimensionless moderately dense (GHD) predictions as a function of coefficient of restitution 
(c) dimensionless dilute (GD) predictions as a function of coefficient of restitution. See legends presented in Figure 1. The 
dimensionless inputs are as follows: mi/m.2 = 8, cri/cra = 2, and 4>i/4> = 0.5. 



coefficient of restitution. In an clastic, moderately dense system (<f> ~ 0.4), results (Fig. 9a) indicate that the prediction 
of thermal conductivity from GHD theory is 4 times larger than that of its dilute (GD) counterpart. For systems of 
lower densities (<fi — 0.1), the discrepancies range from 1% (a = 0.9) to 6% (a = 0.5) (Fig. 9a). 

Similar to the mutual diffusion coefficient, the dilute form of the Dufour coefficient takes on a zero value at certain 
a, thereby making the dense-to-dilute value diverge at this value of a. Because this value occurs at a practical value 
of a = 0.63 (whereas -Dn, dilute diverges at a = 0.52), the dense and dilute predictions of the dimensionless Dufour 
coefficient D* 1 were instead plotted separately against the coefficient of restitution (Figs. 10b and 10c), with the non- 
dimensionalization denned in Eq. (|2.4[) . As expected, the dilute prediction of the Dufour coefficient is independent of 
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FIG. 11: Dufour coefficient {D q ^)- ratio of moderately dense (GHD) to dilute (GD) predictions as a function of (a) overall 
volume fraction and (b) coefficient of restitution. See legends presented in Figure 1. 

the volume fraction (Fig. 10c). 

The differences in magnitude between dilute and moderately dense predictions are non-trivial for both D q \ and 
D q 2- More specifically, the discrepancies that exist between the predictions of D q> \ and Dq^diiute are up to 2 orders 
of magnitude in some cases (Figure 10). For a moderately dense, inelastic system (<j> = 0.5, a — 0.7), the dense-phase 
prediction is over 100 times greater than its dilute counterpart (Fig. 10a). Even at a much lower volume fraction 
(4> = 0.01), the discrepancy between dense and dilute predictions is at least 30%. The differences between dense and 
dilute predictions of D q ^, shown in Figure 11, are less pronounced than D q< x, however, still quite significant. In fact, 
results indicate at least a 20% discrepancy between predictions at a volume fraction (j> = 0.1 (Fig. 11a). As for the 
mutual and thermal diffusion coefficients, the dilute (GD) Dufour coefficient presented in Ref. [25[ is defined using 
different spatial gradients than those used in the dense (GHD) theory and shown in Eq. (|2.2[) . As done before, a 
conversion has been applied to compare the Dufour coefficients by using the same representation for the heat flux. 

V. SUMMARY 

To date, the understanding of particle segregation within polydisperse, rapid granular flows is somewhat limited 
due to a wide array of complexities that arise during the associated derivation of continuum theories. As previously 
mentioned, the two most common simplifications used in previous theories have been a Maxwellian velocity distribution 
and an equipartition of energy. This study focuses on two particular theories, neither of which assumes the above 
conditions. The first was proposed by Garzo and Dufty p5l l2rjj for binary, dilute mixtures (referred to as GD 
theory), and the second was recently proposed by Garzo, Hrenya and Dufty [H, [29[ for binary, moderately dense 
mixtures (referred to as GHD theory) . In order to gauge the importance of this dense-phase extension, the transport 
coefficients and equations of state predicted by GHD theory were compared to their dilute counterparts (GD theory). 
Furthermore, although not the focus of this study, it is worthwhile to mention that the CPU time required to evaluate 
the dense-phase coefficients was typically three times the requirement for its dilute counterpart. 

A systematic comparison was carried out for three different cases (equal size and different mass, equal mass and 
different size, and different size and mass) over a range of mixture parameters (diameter ratio, mass ratio, and volume 
fraction ratio), the details of which are listed in Table 2. Though this study focuses on a case of different-sized 
species with the same material density, similar trends were observed for all other cases analyzed. Results indicate 
that transport coefficients and equations of state predicted by GHD theory are substantially different than those 
predicted by dilute (GD) theory. Also, significant differences between predictions were reported for fairly dilute 
systems (</> = 0.1). In particular, the discrepancy between predictions was found to be as large as an order of 
magnitude. Certain coefficients, namely the mutual diffusion coefficient -Dn, revealed that the magnitude and sign 
were different for the two theories. Naturally, the level of desired accuracy may vary between users of the theories. If, 
for example, 5% deviation between the GHD and GD predictions is deemed acceptable, then the need for a dense-phase 
correction is quite evident since the vast majority of quantities predicted by GHD theory are either larger or smaller 
than GD theory predictions by the 5% limit. Nonetheless, it is worthwhile to point that the comparison of dense- and 
dilute-phase predictions for a binary mixture presented here are independent of flow geometry. It is expected that 
is some flow geometries, one or more of the transport coefficients may dominate, while in other geometries another 
coefficient (s) may dominate. Such differences are system-dependent and should be taken into account when using the 
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results contained herein. 

Given the importance of the dense-phase corrections on the equations of state and transport coefficients, several 
follow-on studies are warranted: application of the theory to segregating systems in order to better understand the 
dominant segregation mechanisms (some previous studies have been carried in the tracer limit 

EMI), 

comparison 

with experimental and/or molecular-dynamics simulation data for purposes of validation, and application of the 
theory to a continuous particle size distribution. It is worthwhile to note that the GHD theory has been incorporated 
recently into the open-source, public MFIX code (https:/ /mfix. nctl.doe.gov/) for the case of binary mixtures, thereby 
increasing its availability to a wider class of researchers. 
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Appendix A: Corrections to previous results 

In this Appendix we explicitly state some changes we have made in the original papers [28| and [1^] to correct several 
errors and/or misprints we have found while working the present manuscript. With these changes, the interested reader 
can easily obtain the complete set of equations for the mass, heat and momentum fluxes and the cooling rate displayed 
along Section |TTJ 

Now, we list the changes affecting both papers: 

• In Eq. (6.18) of Ref. [28| . the term ri£d nj lnx^ appearing in the second line of the right hand side of this 
equation must be replaced by rijd n . \nxu ■ This change affects to Eq. (C6) of Ref. [29| so that, its left hand 



side should read 



(0) d ( 91n Xi° , t i 
i=l \ j / 

Equations (|2 . 7[) and (|2.8j) can be derived after considering these changes. 



A factor "3" and the diameter are missing in the expression (F25) of Ref. [28[ in the collisional contribution 
to the heat flux. Thus, the third and fourth lines of the right hand side of this equation become 



P =i 



where fig = mi/(mi + mj). These changes also affect to Eqs. (7.14)-(7.16) of Ref. [28[ and to Eqs. (3.37)-(3.39) 
of Ref. [29| (collisional contributions to the heat flux transport coefficients). Taking into these changes, one gets 
Eqs. ([2~T7| and (|2~2TJ)) of the present paper. 

The first line of the right hand side of Eq. (3.61) of Ref. [29| must be corrected. It should be given by 

- d + 2 n % n 3 T? ( m, x % u u - (^6g _JLl^^!lx 1 9\nT t \ 

2 m 4 T 2 \pT t ^ me n e T e t] d + 2 n.Tf dn, 1 ,IT, On, ) ' ' " 

Equation (|B10[) of the present paper can be obtained after this change. 
• In Eq. (2.16) of Ref. [29j . the right hand side should read 

2 

d{d + 2) ni euD 
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A minus sign lacks on the second line of Eq. (2.19) of Ref. [29| . Thus, this line should read 



\- (0) 



.2; 3 -=i 

• In Eq. (2.20) of Ref. (29|, the term e^u should be inside the summation sign. Moreover, the partial densities 
TiiTij must be also included inside the summation sign. Thus, the first line of Eq. (2.20) should read 

c ~ 4dr(^) nT^^^ 6 '' 5 ^ x ^ 

V 2 / 2=1 J = l 

Equation (|2.42j) of the present paper can be easily obtained after taking account these changes for the cooling 
rate. 

• The summation $^=1 is missing on the right hand side of Eq. (3.34) of Ref. (29j . 

• On the right hand side of Eqs. (A28) and (A29) of Ref. (29|, the term 1 + ctij should be changed to 1 + eta- 

• The ratio me/mi on the left hand side of Eq. (3.21) of Ref. [29[ should be removed. 

• In Eq. (A12) of Ref. J2jj, the factor d + 5 near the end of the second line should be replaced by the factor d + 3. 

• The right hand side of Eq. (3.54) of Ref. HH should read 

A fc = E *? = £> + ••■ 

i=l i=l 

Appendix B: Some explicit expressions 

The kinetic part of the transport coefficients D q i and A are given by Eqs. (|2.15[) and (|2.16[) , respectively. The 
(dimensionless) Sonine coefficients A* are defined by the matrix equation 

( "721 722 ~K* ) ' ( A* ) = ( X ) ' (B1) 



where 



A* = ± - ^^Df + ^ (i) r^M^Ma,)^^) . (B2) 



3=1 



The expressions of the (reduced) collision frequencies 7*- and uj*- can be found in the Appendix A of Ref. [291 ] . 
Moreover, in Eq. (A2) we have introduced the quantity 

A i:j = (d + 2)(M 2 . - 1) + (2d - 5 - 9n, ,;•.!/,,. ,\/,, + (d - 1 + 3a y + (><.;, i.\/f, + 1 + a„-) 2 . (B3) 

The solution to Eq. (Al) is elementary and gives 

A * = (722 - 20^1-712^ \* = (7i*i-2C)A;-72iA; 

1 4(* 2 - 2(7r i +7| 2 )C- 7M +71*172*2' 2 4C* 2 - 2(7^ + 72 * 2 )C* - 7i* 2 7 2 *i + 7i*i7 2 *2 ' 1 ' 

With these results the kinetic part A fe * can be written as 

vfe * _ X(lh - 2C* - 72*1 ) + A 2 *(7i*i - 2C* - 7fa) , / Ti 72 \ nT * 



\/c* "iv la I iii 1 - -i \ 111 / 1^/ 1 / i-i 1 n-i * CB r ) 

4C* 2 - 2( 7 f! + 7 2 * 2 )C* " 7i* 2 72i + 7i*i7 2 * 2 \M 12 M 21 ' ' 
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The kinetic part of the transport coefficients is given in terms of the Soninc coefficients d* By using matrix 



notation, the coupled set of four equations for the coefficients 

l^g,ll! ^<3,12> ^q,2lJ ^9,22} 

can be written as 



a y — v 



Here, is the column matrix defined by the set (A6) and A^,/ is the square matrix 



A, 





721 







7n - | C* 


721 



7l2 



722 " f C* 





7l*2 





72*2 - f C* / 



The column matrix Y is 



Y = 



l q.l2 

T 

q,21 
V ^,22 ) 



where 



mi + 7712 



'Xi^fillj 



d-fj 
drij 



mi+m 2 9 \ ~> 
e=i 



D 



mi 



n 3 c>C (0) 
vq drij 



T d/2 



mi + m 2 



d(d + 2)r(f) m 
1 



Sit + 



d-fi 



1=1 



mi'Yi drij 



-A, 



In Eq. (A10), A i3 - is defined by Eq. (|B3|) and is given by 



B 



m 2 T 2 

- (d + 8)M?- + (7 + 2d - 9ay)My + (2 + d + Sa% - :\a,,)M 2 , + 3M? (1 + ) 2 — ^ 



+ [(d + 2)M 2 + {2d - 5 - !)r»,,i.\/,,.\/„ + (d - 1 + 3a y + Qa^M^] 



-(d + 2) 1 + 



rrijTi 



The solution to Eq. ( A7) provides the expressions of d* y . The result is 



"0,11 — 



4 ^g,2i7i2 - H,n722 + 6d g ,iiC* 
47i*272i + (27i*i " 3C*)(3C* - 2 72 * 2 ) ' 



9,12 



H-22712 ~ 4rf 9,1272*2 + 6dq, 12 C* 

47i*27 2 *i + (27ii " 3C*)(3C* - 2 72 * 2 ) ' 



(B6) 
(B7) 

(B8) 



(B9) 



(BIO) 



(Bll) 



(B12) 



The kinetic part D*\ can be easily obtained when one takes into account Eqs. (|2.15|) and (A10). The result is 



k * _ 471*2(^,21 + d,,a 2 ) + 2(d,, u + d 9 ,i 2 )(3C* " 272*2) , / 7i 



47i*272*i + (27i*i-3C*)(3C- 2 72 * 2 ) 



VMi 



72 

M 2 i 



(B13) 



The expressions of d* 22 , d* 21 and D^* 2 can be obtained from Eqs. (|B12[) and (|B13|) . respectively, by interchanging 



l-H-2. 



In order to get the dependence of the transport coefficients on the parameters of the system, one needs to know the 
explicit forms of Xij an d fJ-i- For hard disks (d = 2), a good approximation for the pair correlation function Xij is [13] 



Xij 



16(l-0) 2 a l3 M 2 



(B14) 
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where 



M„ = J2 x s a r s 



(B15) 



s=l 



The expression of the chemical potential \ii of the species i consistent with the approximation (IB14[) is [3£ 



T 



ln(A 2 n 4 ) - ln(l - 0) + 



4Mo 



90 



0(1 - 100) 8 



+ ln(l - 0) 



Mf (1 - 0) 2 M 2 1 - M; 



ln(l - 0) 



(B16) 



where Aj(T) is the (constant) de Broglie's thermal wavelength |37l|. In the case of hard spheres (d = 3), we take for 
the pair correlation function Xij the following approximation |38j 



1 



Xij 



<Ti<TjM-2 1 



1-0 2(l-0) 2 CTijMs 

The chemical potential consistent with (|B17| is [13] 



2 (1- 



<J lJ M 3 
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Appendix C: Expressions for a low-density granular binary mixture 

In this Appendix we include the explicit expressions of the transport coefficients and the cooling rate for a dilute 
binary granular mixture. These expressions can be easily obtained from the results derived in Scc.|lI]for a moderately 
dense binary mixture by taking the limit na 2 — > 0. 



1. Mass and heat flux transport coefficients 

The expressions of the reduced coefficients Df* , , and D 12 are given by 
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where is given by Eq. (|2.9p with X12 = 1- Upon deriving Eqs. (|C1|) - (|C3|) . use has been made of the identities 
p* = 1 and dp/dn,i = T. 

The collisional transfer contributions to the heat flux transport coefficients vanish in the low density limit (na 2 — > 0). 
Thus, only their kinetic contributions must be considered. In dimcnsionless forms, the thermal conductivity A* and 



the Dufour D* { coefficients are 
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^.i = di n + < 21 + - ^) ^l-Dii) D*, 2 = d; 22 + < 12 + - ^ z^, (C5) 

where the coefhcients A* and d* ^ are given by Eqs. (|B4j) and (|B12|) . respectively, with 

T * mi+m 2 2^ fx u *i ~ £* Si i r,*r\ 

= ^— X ^ g [ Sij ~ , rj ) ' (C6) 

-* _ mi + m 2 gTj mi + m 2 2 A q£ - Cfa rc,- dC (0) x . 

The dilute forms of the collision frequencies oj*j and 7*j can be obtained from their dense counterparts (Appendix A 
of Ref. [29[) by simply taking Xij = 1- 

2. Pressure tensor 

In the low-density limit, the hydrostatic pressure p = nT . the bulk viscosity k = and the shear viscosity rj has 
only kinetic contributions. It is given by rj = {p/vo)v* where rj* = r]* + ?? 2 . The partial contributions 77* are given by 
Eqs. (|2~30l) and (|23T|) with 

?h= 2:171 > '? 2 = 2; 2 7 2 . (C8) 

As before, the (reduced) collision frequencies T*j can be easily obtained from their corresponding dense forms by 
considering \ij = 1- 

3. Cooling rate 

The first-order contribution ( u to the cooling rate C vanishes in the dilute case [see Eqs. (|2.40p - (|2.45p ]. The dilute 
expression for the coefficient f^ * 1 = can be obtained from Eq. (|2.38p by taking Xij = 1- 



[1] S. Savage, in Developments in Engineering Mechanics, edited by A. P. S. Selvadurai, Elsevier, Amsterdam, 1987, pp. 
347-363. 

[2] J. Bridgwater, in Granular Matter: An Interdisciplinary Approach, edited by A. Metha, Springer, Berlin, 1994, pp. 161-194. 

[3] J. Ottino, D. Khakhar, Mixing and segregation of granular materials, Ann. Rev. Fluid Mech. 32 (2000) 55-91. 

[4] IUTAM Symposium on Segregation in Granular flows, edited by A.D. Rosato, D.L. Blackmore, Kluwer, Dordrecht, 2000. 

[5] C. M. Hrenya, Kinetic theory for granular materials: Polydispersity, in Computational Gas-Solids Flows and Reacting 
Systems: Theory, Methods and Practice, edited by S. Pannala, M. Syamlal, T. OrBrien, IGI Global, Hershey, PA. In press. 

[6] S. Hsiau, M. Hunt, Granular thermal diffusion in flows of binary-sized mixtures, Acta Mechanica 114 (1996) 121-137. 

[7] J. J. Brey, M. J. Ruiz-Montero, F. Moreno, Energy partition and segregation for an intruder in a vibrated granular system 
under gravity, Phys. Rev. Lett. 95 (2005) 098001. 

[8] J. Galvin, S. Dahl, C. Hrenya, On the role of non-equipartition in the dynamics of rapidly flowing granular mixtures, J. 
Fluid Mech. 528 (2005) 207-232. 

[9] V. Garzo, Segregation in granular binary mixtures: Thermal diffusion, Europhys. Lett. 75 (2006) 521-527. 
[10] M. Schroter, S. Ulrich, J. Kreft, S. B. Swift, H. L. Swinney, Mechanisms in the size segregation of a binary granular 

mixture, Phys. Rev. E 74 (2006) 011307. 
[11] R. Wildman, J. Jenkins, P. Krouskop, J. Talbot, A comparison of the predictions of a simple kinetic theory with experi- 
mental and numerical results for a vibrated granular bed consisting of nearly elastic particles of two sizes, Phys. Fluids 18 
(2006) 073301. 

[12] D. K. Yoon, J. T. Jenkins, The influence of different species granular temperatures on segregation in a binary mixture of 

dissipative grains, Phys. Fluids 18 (2006) 073303. 
[13] X. Liu, M. Metzger, B. Glasser, Couette flow with a bidisperse particle mixture, Phys. Fluids 19 (2007) 073301. 
[14] V. Garzo, Brazil-nut effect versus reverse Brazil-nut effect in a moderately granular dense gas, Phys. Rev. E 78 (2008) 

020301(R). 

[15] V. Garzo, Segregation by thermal diffusion in moderately dense granular mixtures, Eur. Phys. J. E 29 (2009) 261-274. 



21 



[16] V. Garzo, F. Vega Reyes, Mass transport of impurities in a moderately dense granular gas, Phys. Rev. E 79 (2009) 041303. 
[17] J. Jenkins, F. Mancini, Balance laws and constitutive relations for plane flows of a dense, binary mixture of smooth, nearly 

elastic, circular disks, J. Appl. Mech. 54 (1987) 27. 
[18] V. Mathiesen, T. Solberg, H. Arastoopour, B.H. Hjertager, Experimental and computational study of multiphase 

gas/particle flow in a CFB riser, AIChE Journal 45 (1999) 2503-2518 (1999). 
[19] L. Huilin, D. Gidaspow, E. Manger, Kinetic theory of fluidized binary granular mixtures, Phys. Rev. E 64 (2001) 061301. 
[20] M. F. Rahaman, J. Naser, P.J Witt, An unequal granular temperature kinetic theory: Description of granular flow with 

multiple particle classes, Powder Technol. 138 (2003) 82-92. 
[21] H. Iddir, H. Arastoopour, Modeling of multitype particle flow using the kinetic theory approach, AIChE Journal 51 (2005) 

1620-1632. 

[22] J. Jenkins, F. Mancini, Kinetic theory for binary mixtures of smooth, nearly elastic spheres, Phys. Fluids A 1 (1989) 
2050-2057. 

[23] B. Arnarson, J. T. Willits, Thermal diffusion in binary mixtures of smooth, nearly elastic spheres with and without gravity, 

Phys. Fluids 10 (1998) 1324-1328. 
[24] S. Chapman, T. G. Cowling, The Mathematical Theory of Nonuniform Gases, Cambridge University Press, Cambridge, 

1970. 

[25] V. Garzo, J. W. Dufty, Hydrodynamics for a granular mixture at low density, Phys. Fluids 14 (2002) 1476-1490. 
[26] V. Garzo, J. M. Montanero, J. W. Dufty, Mass and heat fluxes for a binary granular mixture at low density, Phys. Fluids 
18 (2006) 083305. 

[27] V. Garzo, J. M. Montanero, Navier-Stokes transport coefficients of d-dimensional granular binary mixtures at low density, 
J. Stat. Phys. 129 (2007) 27. 

[28] V. Garzo, J. W. Dufty, C. M. Hrenya, Enskog theory for polydisperse granular mixtures. I. Navier-Stokes order transport, 
Phys. Rev. E 76 (2007) 031303. 

[29] V. Garzo, C. M. Hrenya, J. W. Dufty, Enskog theory for polydisperse granular mixtures. II. Sonine polynomial approxi- 
mation, Phys. Rev. E 76 (2007) 031304. 

[30] A. Goldshtein, M. Shapiro, Mechanics of collisional motion of granular materials. Part 1. General hydrodynamic equations, 
J. Fluid Mech. 282 (1995) 75-114. 

[31] J. J. Brey, J. W. Dufty, A. Santos, Dissipative dynamics for hard spheres, J. Stat. Phys. 87 (1997) 1051-1066. 

[32] M. Lopez de Haro, E.G.D. Cohen, J. Kincaid, The Enskog theory for multicomponent mixtures. I. Linear transport theory, 
J. Chem. Phys. 78 (1983) 2746-2759. 

[33] H. van Beijeren, M. H. Ernst, The non-linear Enskog-Boltzmann equation, Phys. Lett. A 43 (1973) 367-368; The modified 
Enskog equation for mixtures, Physica 70 (1973) 225-242. 

[34] V. Garzo, J. M. Montanero, Shear viscosity for a moderately dense granular binary mixture, Phys. Rev. E 68 (2003) 
041302. 

[35] V. Garzo, J. W. Dufty, Homogeneous cooling state for a granular mixture, Phys. Rev. E 60 (1999) 5706-5713. 
[36] A. Santos, private communication. 

[37] T. M. Reed, K. E. Gubbins, Applied Statistical Mechanics, McGraw-Hill, New York, 1973, Chap. 6. 

[38] T. Boublik, Hard-Sphere equation of state, J. Chem. Phys. 53 (1970) 471-472 (1970); E. W. Grundke, D. Henderson, 
Distribution functions of multi-component fluid mixtures of hard spheres, Mol. Phys. 24 (1972) 269-281; L. L. Lee, D. 
Levesque, Perturbation theory for mixtures of simple liquids, Mol. Phys. 26 (1973) 1351-1370. 



